A consistent formalism for the Thomas-Ehrman Level Displacement 
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Usage of the Thomas-Ehrman Level Displacement formalism has been examined. Mistakes and 
inconsistencies are found in several papers, being repeated in subsequent works. Here, we present a 
complete formalism with a consistent set of definitions. Full algorithms are made available, both as 
a FORTRAN source file and as a user-friendly Visual Basic executable tool, available for download on 
the World Wide Web. 
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I. INTRODUCTION 

The Thomas-Ehrman Level Displacement formalism 
(TELD) 0,0 is an established technique for calculat- 
ing the level displacement between mirror pairs. It is 
found to be particularly useful in situations where a re- 
action proceeds via a proton resonant state in a proton- 
rich nucleus. Largely, this usefulness derives from the 
fact that such states are above the particle decay thresh- 
old, usually resulting in proton partial widths too narrow 
to be measured experimentally. Thus, by appealling to 
the charge symmetry of the nuclear force, one may make 
use of relatively abundant spectroscopic data of analogue 
states in the mirror nucleus to determine the properties 
of the astrophysically interesting states. Examples in 
the literature are the ^iNe-^^Na Q, ^op-^ONa [3, ^^O- 
i^Ne 0, 22Ne-22Mg ^ and ^STi-^eCr 0,11 mirror nu- 
clear pairs. However, a survey of the literature finds in- 
consistency in the definition of critical parameters, lead- 
ing to errors in the calculations. In the present work 
a complete and consistent TELD formalism is presented 
and made available for wider use. 



II. THE WAVE FUNCTION 



Here we reproduce and expand upon the original work 
of Thomas [2'|], using, for consistency, exactly the same 
terminology. The channel radius of two interacting bod- 



ref 0) may be written that satisfies the wave equation 
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ies is defined as Oc = 1.44 X {Al^^" + ^2^") fm, with Ale 
and A2c being the mass numbers of the bodies of the 
pair; the reduced mass is Mc = AicA2c/{Aic + A2c); the 
energy of relative motion is Cc, which may be positive or 
negative. The subscript c is used to describe all of the fea- 
tures of the channel, unless it is necessary to distinguish 
the positive-energy (ec+ > 0) from the negative-energy 
(ec- < 0) channels in which case the symbols c+ and c— 
are used, respectively. 

For external wave functions, a radial factor (Equ. 1 of 



(1) 



where a prime signifies differentiation with respect to r 
(in the following descriptions, all the derivatives are with 
respect to r unless stated otherwise). The interaction 
potential may be written 



Oc = ZicZocS-r;. 
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where the nuclear potential term disappears in the ex- 
ternal region. In the notation of Yost, Wheeler, and 
Breit the positive-energy solution, which is regular 
at the origin, is designated by F{kr) and has the asymp- 
totic form for large r, 

Fc+ ~ sin(a; — — i?7r — ?7ln2a; + a). (3) 

Likewise, there is a solution which is linearly independent 
of F and irregular at the origin which is conveniently 
taken with the asymptotic form for large r, 

Gc+ ~ cos(a; — ^^tt — rj\ii2x + cr). (4) 
The quantities entering Equ. [3] and 3] are 

Xc± — kr 

kc± = p/h={2M,\e\/hY^^, 
with Sommerfeld parameter 

r]c± = McZicZ2ce^ /h^k ZicZ2ce^ /hv, 

and 

(Tc+ = argr(l + ^ + ^77). 

It is worth noting that x is replaced with p in some for- 
mulations. 

The general solution of this equation, F{r), is a linear 
combination of F and G. The Wronskian relation for 
these two particular solutions, which directly follows from 
Equ. [1]- His 
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c+- 



(5) 



2 



Extensive tables [l0| and several computer codes [ll|, [T^, 
[Tsj have been developed for evaluating F and G and their 
derivatives when 77 > 0. 

I 



Whittaker function and its derivative may be accurately 
calculated using the whittaker_w 15] computer code. 
However, it is useful to note that if there is no Coulomb 
interaction in a c— channel, one has from Equ. [S] for s, 
p, d, and / orbitals the simpler relations 

W^^.{2x) = e-^ (7) 
W^^.{2x) = (l + x-')e-'^ (8) 
Wq^'^{2x) = {l + ix-^ + 3a;~2)g-^ (9) 
Wq^7_{2x) = (1 + Qx-^ + 15a;"2 _^ i5a;-3)e-^ (10) 

which can be used for checking the results from a more 
complicated code. 

In discussing conditions at the nuclear surface, one 
needs to evaluate the real and imaginary parts of the 
logarithmic derivatives, gc = E' /E, and these are 0], 



g^^ = {FF' + GG'){F^ + G^)-! (11) 

g^^™ = fc(^^2^G2)-i (12) 

g^^ = W'W-^ (13) 

5c'" = (14) 



where gc — g^'^ + i^^™ and rr — Qc- Although the sim- 
ple WKB approximation [3, 0] can perform well in cal- 
culating the logarithmic derivatives of the Coulomb and 
Whittaker functions in specified regions, modern com- 
puter codes perform essentially exact calculations and 
are preferred. For example, the difference between the 
WKB approximation and the exact evaluation of gc-, 
performed using the code whittaker_w 15], in the region 
0.1< X <5.0, 0.1< 77 <10.0 (roughly corresponding to 
■^ic-^2c <20 and bound nucleon energy 0.1< \Ei,\ <10.0 
MeV), can be as large as 3% (£=1), 1% {£^2) and 0.6% 
(£=3), respectively, and becomes smaller as ry increases. 
In the case that £=0, the differences can be as much as 
10% when x <0.1 and 77 <0.9, or X <0.2 and rj <0.5, 
or X <0.5 and 77 <0.3. In addition, it should be noted 
that a frequently used subroutine, COULFG |l2|, was un- 
able to reproduce the results (F, G and their derivatives) 
of the subroutine RCWFN near values of 7; = 1.70a; -I- 5.13 
unless the ACCUR variable was set to be less than 10~^^ 
(of course the smaller the safer, e.g., 10"'^'' if possible). 



For the c— channels, only the solution to Equ. [1] van- 
ishing at large distances from the origin, can occur; it is 
the Whittaker function [13, [II, 



(6) 
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III. CALCULATION OF LEVEL 
DISPLACEMENTS 

Under the assumption of charge symmetry of nuclear 
forces, the nn and pp nuclear interactions are identical. 
The difference in the excitation energy of levels in mirror 
nuclei is therefore due to differences in the Coulomb, elec- 
tromagnetic spin-orbit, and mass energies. As suggested 
by [3], one can evaluate this by considering Hp — Hn = V. 
Furthermore, in the one-level approximation, irrespective 
of the boundary conditions, in the internal region the pro- 
ton and neutron wavefunctions are the same to within a 
multiplicative constant. Under these assumptions, one 
obtains (Equ. 30a in ref. 

E„-E.p = -{V)r + Axn-Axp, (15) 

where there is a boundary condition that satisfies 6„c — 
bpc (defined in Equ. 24c of [3]). Here, {V)r is the mean 
value of V in the internal region, and En and Ep are 
the Eigenvalues satisfying i?„(p)*„(p) = £^„(p)*„(p) of 
the nucleus with the odd neutron or proton, respectively. 
The difference of the level displacements of the neutron 
and proton states (see Equ. 30b of [7]), 

Aa = Axn -Axp^~J2 ilig^c - gpcl, (16) 

is referred to as the boundary condition level displace- 
ment. The energy difference of corresponding levels of 
mirror nuclei, following from Equ. [15] and 1161 can be 
written as 

{E* E^-) - {e; E^p^-) = Al- A{- , (17) 

where it is assumed that the quantity {V)r is the same 
in the excited state as in the ground state. The scripts 
* and g.s. denote the corresponding quantities are be- 
ing evaluated with respect to the excited state and the 
ground state, respectively. Assuming that the level dis- 
placements of the two ground states are the same simpli- 
fies this relation further, 

E*{n)-E*{p)^Al, (18) 

where E* (n) and E* (p) are the corresponding excitation 
energies in the mirror nuclei. Therefore, the observed 
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energy difference between mirror nuclear states is due to 
different level displacements, A;^„ and A\p, in the two 
nuclei. 

We find that the definitions used by various authors 
of partial width Txc, reduced width j'^^, dimensionless 
reduced width 9^ and dimensionless single-particle re- 
duced width 6l2p of a level A, are not consis- 
tent, as shown in Table [11 F ollowing the previous work, 
the definitions of French [l0| are adopted in the present 
work, though it is possible to achieve consistency using 
the alternative definitions 0, ll^L,One should pay atten- 
tion that the definition of 7^^ 19] is different from that 
of by a factor of S/oc- Thus, for the positive-energy 
channel. 



S„=13.1f 



Eb=6.765 



E*(n) =6.424 



E*(p)=6.246 



0.0 



jE,. 


= 1.356 




=4.890 




0.0 



46n 



Ti(^^Ti+ii) 



46/ 



Cr(*'V+p) 



(19) FIG. 1: Definitions and calculation results of relevant quanti- 
ties in case of the ""'Ti-^^Cr mirror pair. The calculated values 



with the Coulomb penetrability Pc — x/{F'^ + G"^); and 
for the negative-energy channel, 



. 0^ (20) 



Where F^,G'^ and represent differentiation with re- 
spect to X. Equ. [m & [20] were defined as and 
in previous literature by assuming 6„c = bpc = (see 
Fig. [ID. 

The reasonable assumption that the reduced widths, 
7^, are the same for the mirror levels leads to an as- 
sumption of 9l = 9p = 9"^. Thus, the excitation-energy 
displacements of mirror nuclei can be expressed as 



(in units of MeV) are with respect to the analogue states at 
6.424 MeV (in ^'^Ti) and 6.246 MeV (in ""^Cr) [Sj. 



The final ingredient which is needed to allow calcu- 
lation of the level shift is the dimensionless reduced 
width, 9^. This can be calculated according to the 
relation 9^=0^39'^ [H, where the dimensionless 

single-particle reduced width 9'^p has already been de- 
termined [HjHSIj ^'^d the factor C^S is calculated via a 
shell-model code such as OXBASH [H. 



^ 2Mcal 



9l [P,{FF', + GG^: 



' W 



\E=Er- 

\ " / \E=Ei 
(21) 

where Er [=E*{p) - S^,] and E^ [=\E*{n) - S'^^W are the 
energies relative to the respective nucleon thresholds {S^ 
and 5^; the superscripts p and n refer to the odd-proton 
(or proton-rich) and odd-neutron (or neutron-rich) nu- 
clei, respectively; the subscripts p and n denote the cor- 
responding nucleon separation energies). A pictorial rep- 
resentation of the physical meanings of the parameters 
described here are illustrated in Fig. [T] which shows the 
case of the 6.424 MeV state in ^^Ti and its analogue state 
in '*^Cr 21]. It is clear from the figure that the level dis- 
placement may equally well be written as, 

E*{n) - E*{p) = [E*{n) - S^] - [E*{p) - S'^] 



- Er, 

(22) 

where the quantity E"^ is different from that of Eb defined 
above. It appears that confusion over these definitions 
has, in part, been the source of errors in the past. For ex- 
ample, in several previous works [4, 5, 6,'7j the level dis- 
placement has been written as E*{n) — E*{p) = Eb — Er 
in contrast to the correct expression of E*{n) — E*{p) — 
E'j^ — Er which follows from Marion et al. [3]. However, 
the pressent work validates that the correct expression 
was used in the calculations (221] . 



IV. ONLINE RESOURCES 



The algorithms developed here have been made avail- 
able online [2^. The Visual Basic tool requires the 
TELD_VB and TELD executables be in the same directory. 
Opening the TELD application provides a window in which 
one defines the parameters of the mirror states to be con- 
sidered [A, Z, 5„, Sp and E.j.{n)). One also could change 
values for the channel radius, the orbital angular momen- 
tum of the single particle and spectroscopic factor for 
the state. Activation of the "calc" button then performs 
the TELD formalism and returns the resulting excitation 
of the proton-rich analogue state and the proton decay 
width for this state. This last parameter is calculated 
using the relation 



ih^PcOl 



(23) 



In addition, a Fortran source code, which allows users to 
make changes whenever required, is also encloded in the 
TELD.rar package. Figure [2] shows a screenshot from the 
TELD program with all relevant parameters for the case 
of the 6.424 MeV state in ''^Ti and its analogue state in 
46 Cr t21|]. 
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TABLE L Various definitions of TactL, and 6% [ig 



of a level A in the literature. 

fl2 



Thomas [2] 
Lane and Thomas [18. 20] 
French [19] 



2Pc(7L/ac) 



2Maac 



^ (r)dr 



li-^ {r)dr 



2 (r)»iT 



lELD vl.O 



(e. g. , 45Ti+ii) 
Znio AnIT 



Zl |22 M |45 
SnOfleV) |13. 189 



Other IiLfomati on- 



. r- 

rO(£m): )l.44 
Option: [i 



e 2 

Exh): |6. 424 (WeV) 



(e. g., 45V+p) 
Zp [i Ap [T 



Z2 |23 A2 |45 
SpftleV) |4.890 



Output - 



EkCp): |6.246 (WeV) 
|0.218E+00 keV 
r^* |0.218E+01 keV 
Aj |0.4e6E+00 keV 
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SUMMARY 



A complete and consistent Thomas- Elirman Level Dis- 
placement formalism has been presented and made avail- 
able on the World Wide Web. With this, if one has 
knowledge (or makes a reasonable assumption) of the 
quantity 9^ and spectroscopic factor 5", one may esti- 
mate the location of the mirror to a known excited state. 
Alternatively, experimental measurement of the mirror 
level displacement provides a route to determining 6^ (or 
S factor). 
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